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We give a technique for calculating the occupation number of quantum fields in time-dependent 
backgrounds by using the relation between one-dimensional quantum oscillators and two-dimensional 
classical oscillators. We illustrate our method by giving closed analytical results for the time- 
dependent spectrum of occupation numbers during gravitational collapse in any number of dimen- 
sions. 



Frequently we are interested in calculating the spectrum of excitations of a quantum field produced due to a time- 
dependent background. The context can range from nuclear physics, condensed matter systems, gravitational systems, 
and cosmology, and includes particle production during gravitational collapse or in an expanding universe. Usually 
such problems are solved by the method of Bogoliubov transformations [1, 2], and require knowledge of all the modes 
of the quantum field at the initial and final times. Then the excitation spectrum, as well as local information, e.g. 
energy-momentum densities, can be obtained. However, in certain cases it may be difficult to find all the modes of 
the quantum field, especially in a time-dependent situation, and it is desirable to find a short cut to the spectrum of 
excitations, even if all the detailed local information is not available. In this paper, we propose such a method based 
on the functional Schrodinger equation that can be applied to certain systems. We illustrate the method in the special 
case of quantum excitations produced during gravitational collapse. 

Consider a scalar field x) that is free but for interaction with a time-dependent background e.g. a time-varying 
mass parameter or background metric. We can always expand the field in a complete basis of functions, 

$(t,x) = 5^a fc (t)/ fc (x) (1) 

k 

where k denotes a complete set of mode numbers. The f k s need not be a set of (instantaneous) mode functions; the 
only requirement is that any function of x should be expressible as a linear sum of fk's. On inserting this form in the 
action, we obtain 



S= dt 



^a k M k i(t)ai - ]^a k V k i(t)ai 



(2) 



where M and V are real, symmetric matrices that are given by spatial integrals over products of the basis functions 
and their derivatives. They depend on time due to the interaction of $ and the background, and a sum over repeated 
indices in Eq. (2) is implicit. In general, the elements of M will have complicated time dependence and the problem 
appears intractable. However, progress can be made if the time dependence can be pulled out, 

U = W) M (3) 

where M is independent of time. Such a simplification occurs for gravitational collapse because B — > and then M 
has the form in Eq. (3) to leading order in 1/B [5]. Similar simplification will also occur in cosmological spacetimes 
in which all the time dependence is in an overall scale factor. 
With Eq. (3) one can write, 



drj 



-a' k M kl ai - ^rj(^ a k v ki(v)ai 



(4) 



where primes denote derivatives with respect to a new time coordinate r\ that is defined by 

Now we attempt to do a "principal axis transformation". As is standard procedure [3], first we diagonalize M. Then 
by a rescaling transformation, M can be transformed into the identity matrix. Since M is taken to be independent of 
time, the diagonalization and rescaling transformations are time independent.. The transformed potential matrix will 
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still be time dependent. The tractable case is when the transformed potential matrix can be diagonalized by a time 
independent transformation. The final action takes the form 



V 

2 k 



Ml) ,2 

2B( V ) k 



(6) 



where bk 's are the transformed mode coefficients and Vk (jl) are the eigenvalues of the transformed potential matrix. 
It should be noted that the principal axis transformation can be implemented even if the matrices M and V do not 
commute. The theorem that only commuting matrices can be simultaneously diagonalized applies when we restrict 
the diagonalization procedure to similarity transformations. The principle axis transformation involves a rescaling 
which is not implemented by a similarity transformation. 

We can now quantize each mode separately. The Schrodinger equation for a mode is 



1 ° 2 1 2, xrf 



d 



(7) 



where we have omitted the mode label k for convenience and written the equation in a form reminiscent of the simple 
harmonic oscillator (SHO). The only time dependence is due to a time varying frequency uj(rf). 
The Schrodinger equation has the "ground state" solution [4] 
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(8) 



where p^ denotes derivative of p(rj) with respect to 77, and p is given by the real solution of the ordinary differential 
equation 



p., n +uj 2 (t])p = 



1 



(9) 



Initial conditions for p are chosen so that the wave function yields the ground state of the time-independent SHO at 

V = 0, 



P(0) 



1 



where cjq = w(0). The phase a is given by 



a(r)) 



Pv(0) = 



dr]' 



2 Jo P 2 (V') 



(10) 



(11) 



The wave function in Eq. (8) can now be expanded in the "instantaneous" SHO eigenfunctions with frequency oj{rj). 
The occupation number of the mode is given by the expectation value of the quantum number of the excitation. The 
occupation number as a function of time and frequency can be written as [5] 



N(t,LJ ) 



1 



4wp" 



[(-P 2 



1) 



(PPv 



(12) 



where the ui on the right-hand side refers to u>{rf). In (12) we have corrected the pre-factor stated in Ref. [5], 1/4 
instead of 1/^2 [17]. 

The crucial step that remains is to solve Eq. (9) and it is here that the remarkable connection between quantum 
SHO in 1 dimension and classical SHO in 2 dimensions comes into play [6-10]. The solution for p is given by the 
solutions for two classical SHO's £ and x> each of which satisfies the equation 



f vv + u 2 (v)f = 0- 



Then the solution for p is 



UJ 



(13) 



(14) 
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The initial conditions for p can be satisfied provided 

£(0) = 1, 6,(0) =0, X(0) = 0, x v (0)=u o (15) 
The time derivative of p, which occurs in the wave function, can be written in terms of £ and \ an d their derivatives, 

Pr, = — (££tj + XXr))- (16) 

An explicit scheme for constructing solutions with the correct initial conditions can now be described. We start 
with the classical SHO equation (13) and let fa and fa be any linearly independent solutions. Then the solutions with 
the correct initial conditions in Eq. (15) are given by 

t{uo,v) = W~ 1 (fa v (uJo,'n i )h(^o,v) - fi v (^o,Vt)h(^o,v)) (17) 

x(u ,v) = -^oW~ 1 (fa(u} ,rj i )fa(uJo,ri) - fi(va,Vi)M u o>V)) ( 18 ) 

^(cjo,??) = W^ 1 (f 2v (oj , Vi)f iv(^o,v) ~ hv( UJ f),Vt)h v {^o,v)) (19) 

Xt,(wo,v) = -u W~ 1 (fa(u ,ri i )f lr ,(uo,r)) - fa{ua- > r)i)far,{ua,rj)) (20) 
where loq = w(0) and W is the Wronskian defined as 

W = f 2v (u , Vi)f l{uo,Vi) ~ flr,{ u} 0,rii)f2(^0,'ni) ( 21 ) 

Now it is simply a matter of finding p and p v from equations (14) and (16), and then inserting into the occupation 
number formula (12). 

We will illustrate the above method by applying it to particle production of a massless scalar field during gravita- 
tional collapse of a spherical domain wall in 3 spatial dimensions. The function B(t) in Eq. (3) is given by [5] 

B(t) = 1 - 4^ (22) 

where R(t) is the radius of the wall at time t, and Rs is the Schwarzschild radius. For convenience we will set R$ = 1. 
Using the classical solution for R(t) in the late time limit and Eq. (5), a model for w(ry) can be taken to be [18] 

w(f7 ) = = (23) 

Then (13) can be solved in terms of Bessel functions e.g. see [11]. With the help of Eqs. (17-20), we get 

Till 

£ = — [y (2wo)Ji(u) - J (2wb)yi(u)] (24) 
X = — [y 1 (2w )Ji(«) - Ji(2wb)yi(«)] (25) 

= -7rwg[y (2wb)Jb(«) - Jo(2w )^o(w)] (26) 

X, ; = -7rwg[l r 1 (2w )Jo(«) - Ji{2u )Y (u)] (27) 

where u = 2woVl ~~ ^ The particle occupation number per mode at time i and frequency loq is found from Eq. (12) 
and can be written as [5] 



A{e+x 2 ) 



(28) 
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FIG. 1: Occupation number, N, as a function of uiq for rj = 0.92, 0.98, and 0.99 (lowest to highest curve). 



where £, x, Xri are given in Eqs. (24), (25), (26), and (27). 

In Fig. 1 we plot the spectrum of occupation numbers at three different times. The plot can also be made in terms 
of fi, the final frequency corresponding to the original time parameter i; SI is also the frequency relevant for an 
asymptotic observer. The relation between t and r\ is obtained from Eq. (5), 



t = 



drj 



drj 
~1 



= -ln(l-7?). 



and the relation between Q and loq is 



(29) 



(30) 



The occupation number spectrum in terms of fl is shown in Fig. 2, and Fig. 3 illustrates the time dependence of the 
occupation number at three different frequencies. 
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FIG. 2: Occupation number, N, as a function of Q for rj = 0.92, 0.98, and 0.99 (lowest to highest curve at origin). 



Limiting forms of the occupation numbers in different frequency ranges are obtained by using expansions of the 
Bessel functions. The final results can be written as 

(V/ 2 /4, 2ft«e~ t / 2 
N w I 1/(4ttQ), e-*/ 2 <2r!<l (31) 
ll/(64ft 2 ), l«2fl. 
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It is also possible to give a closed formula for the occupation numbers for t — > oo that is valid at all frequencies 

2 f 



N 1 *° 



4(J?+1?) 



(32) 



where all the Bessel functions have argument 20. 

The occupation number in the intermediate regime e~*/ 2 = y/1 — rj -C 2fi -C 1 correspond to the low frequency 
part of the thermal blackbody temperature 

= (33) 

if we take the temperature to be 

in units where the Schwarzschild radius Rg — 1. This temperature is precisely the Hawking temperature [12]. The 
very low frequency limit in Eq. (31) only coincides with the divergent thermal value if we take rj — > 1 (t — > oo). 

The high frequency limit in Eq. (31) does not match the exponentially (Boltzmann) suppressed thermal expression. 
This can be understood by noting that our calculation yields occupation numbers during the collapse and includes 
all excitations that are produced, including transients and non-radiative excitations (vacuum polarization). Such 
non-thermal features are also seen in a treatment of gravitational collapse in a 1+1 dimensional model [13, 14] and 
are separated out by using the explicit mode functions to calculate the energy flux at spatial infinity [19]. Only the 
radiative part, i.e. the Hawking radiation, is thermal and hence exponentially suppressed at high frequencies. If we 
are interested in asymptotic signatures of gravitational collapse, it is indeed only the radiative excitations that are of 
interest. However, if we are interested in the dynamics of the collapse, then it is necessary to consider all modes of 
the scalar field and the full wave function as given in Eq. (8). 

The above comments highlight the advantages and limitations of the technique developed in this paper. For systems 
in any number of dimensions that have the simplifying properties discussed around Eqs. (3) and (6) we can obtain the 
complete wave function for the mode coefficients. The solution for the wave function, Eq. (8) , completely characterizes 
the system and allows us to calculate the occupation numbers analytically and in closed form, even without explicit 
knowledge of the mode functions. However, it becomes necessary to find the mode functions if one requires spatial 
information about the excitations that are produced. Additionally, we have neglected backreaction of the radiation 
on the collapse. This is a difficult problem that has not seen any clear resolution in the literature; some attempts in 
the present formalism that lead to a modification of the wave function describing gravitational collapse may be found 
in [15, 16]. 

We thank Rong Chen, Yi-Zen Chu, Paul Davies, and Dejan Stojkovic for discussions. This work was supported by 
the DOE at ASU. 
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